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The phase diagram of a non-relativistic fermionic system witli imbalanced state populations in- 
i^ ' teracting via a short-range S-wave attractive interaction is analyzed in the mean field approxima- 

^^ ' tion. We determine the energetically favored state for different values of the mismatch between 

'^v , the two Fermi spheres in the weak and strong coupling regime considering both homogeneous and 

|/-\ ■ non- homogeneous superconductive states. We find that the homogeneous superconductive phase 

^v^ ' persists for values of the population imbalance that increase with increasing coupling strength. In 

the strong coupling regime and for large population differences the energetically stable homogeneous 
phase is characterized by one gapless mode. We also find that the inhomogeneous superconductive 
phase characterized by the condensate A(x) ~ A exp (iq ■ x) is energetically favored in a range of 
values of the chemical potential mismatch that shrinks to zero in the strong coupling regime. 
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I. INTRODUCTION 



i-^ , Fermionic superfluidity driven by S-wave short range interaction is a collective phenomenon characterized by for- 

P5 ■ mation of correlated pairs of half-integer spin particles. This effect is relevant in many theoretical and experimental 

O \ fields. Of particular interest are systems of trapped cold atoms of two fermionic species where the interaction between 

I ^1 . the fermions can be varied over a wide range employing a Feshbach resonance [l| . A general feature of such systems is 

that when the interaction between the fermions is increased the systems evolve from a weak-coupled superconductive 

T— I . BCS state of Cooper pairs to a Bose-Einstein condensate (EEC) of tightly-bound diatomic molecules Q, y{. 

K>" ' In a standard BCS superconductor the chemical potentials of the two fermionic species are equal. A small chemical 

O^ I potential difference 5^ cannot disrupt BCS superconductivity. As a matter of fact the superconductive state with 

equal number densities is energetically favored in comparison with a normal state with a fermionic imbalance. On 

the other hand, as pointed out in Q in the weak coupling regime, BCS superconductivity cannot persist for large 

^^ . values of (5//. Indeed there exists an upper limit for (5/x (the so-called Clogston limit), beyond which the homogeneous 

\^ ' superconductive state is no longer energetically favored. 

^> , Syst ems with u neq ual pop ulations of two fermionic species have been widely studied [1, H, B S S H^ IHl H^ 

"^ ■ m [il El El El El, El iM m ill . For values of <5/i larger than the Clogston limit several non trivial phases have 
2 , been proposed. They include the Phase Separation (PS) state Q, the Breached Pair (BP) superfluidity d El Ell, 
the Deformed Fermi sea pairing (DFSP) Q and the non-homogeneous LOFF [g phase (see ^^ for a review). 



"i \ Recent experiments with trapped cold atoms [2J, |25j indicate that both the BCS and the BEC states can sustain 

^ ' large mismatches between the fermionic populations before turning to the normal unpaired state. For a wide range of 

Q ' values of such a mismatch it turns out that the supcrfluid atoms tend to remain in the inner core of the trap, whereas 

Q . the normal atoms in excess are expelled and form a shell surrounding the center of the trap. These observations 

seem to support the PS scenario with phase sep aration [g between superfluid and normal atoms. However, a clear 

signature of the PS phase is still lacking |24ll26l|. 

^^ ■ The LOFF phase has not yet been observed in experiments with trapped cold atoms |2J, |23, although it was 



Jh \ seen in CeCoIns [271, |28j . The portion of parameter space where this phase is favored depends on the form of the 
condensate. In the simplest case, in the weak coupling limit this region is rather narrow and its detection in cold 
atom experiments might be difhcult '24| . However a definite conclusion can be reached only after a detailed study of 
the space dependence of the condensate and the knowledge of its dependence on the coupling strength. 
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The present study represents a preliminary analysis of this problem. Wc determine, in the mean field approximation, 
the phase diagram of an interacting two species fermion system as a function of the strength of the interaction and of 
their chemical potential difference, paying special attention to the LOFF region. We choose a four-fermion coupling to 
mimic the interaction, neglecting all the effects of the trap in the fermionic Hamiltonian |33l |. Due to these limitations 
we do not attempt to make any comparison with experimental data, lea ving this task to future studies. 

Other analyses have been presented already on the same subject [JM UJl iln^. In particular in Ref. [lj| a similar 
study has been performed for the homogeneous case. We differ from these authors for two reasons. First, we also 
include and discuss the LOFF phase. Second, we work at fixed total density but do not fix the difference in population 
between the two species. Under this respect our study may be seen as complementary to Ref. [Q. Due to these 
assumptions, our phase diagram at T = is given in terms of two parameters, 6^ and the coupling constant, similarly 
to the analysis presented in |15l |. 

One of our results is that the homogeneous phase is energetically favored in a range of values of the chemical 
potential difference that increases with increasing coupling strength. Moreover it turns out that the phase transition 
between the homogeneous and the LOFF phase is of the first order (this result was already known in the weak coupling 
regime and we extend it to the intermediate coupling regime) and there exists a critical value of the coupling constant 
where the LOFF window shrinks to a point. For larger values of the coupling the LOFF phase cannot be realized 
and the homogeneous phase has a first order phase transition directly to the normal phase. Further increasing the 
coupling strength and for large separations of the Fermi spheres the excitation spectrum of the energetically stable 
homogeneous phase is characterized by one gapless mode and the transition between the superconductive phase and 
the normal phase becomes of the second order. 

As for the non-homogeneous phase we restrict our analysis to the condensate characterized by a one- wave oscillation 
A(x) ~ A exp (iq • x). Non-homogeneous condensates characterized by a more complicated space dependence, such 
as that obtained by a superposition of several plane waves are expected to produce deeper minima of the free-energy 
F J29j . Therefore our study only provides a lower bound for the gap and the window in 5^ where the LOFF state is 
energetically favored. 

The stability of the homogeneous phase with different Fermi surfaces that we present in this article has been 
already discussed in Ref. [lj| . We will only consider the solutions of our system of equations which are stable and 
that correspond to global minima of the free-energy. We do not examine the stability of the LOFF phase which 
however is expected to be stable [30J . 

Our paper is organized as follows. In Section ^] we study the homogeneous superconductive phase constrained 
by the conservation of the total number of fermions. In Section IIIII we discuss a LOFF phase characterized by a 
single plane- wave and we determine the width of the LOFF window as a function of the coupling. In Section Hvl we 
summarize our results. 

II. MODEL AND HOMOGENEOUS PHASES 

We shall analyze a system comprising fermions of two species, with mass m, interacting via a short range (S-wave) 
four fermion interaction. The Hamiltonian density is 



H= V" ^Jcr{x) [-- /icr ) l/jcrix) - Glpl{x)tp2ix)lp2{x)^pl{x) , 



(1) 



where G > is the four fermion coupling constant. The chemical potentials of the two species can be written as 
fii = fi + Sfi and fj,2 = fJi — Sfi, so that fi is the average of the two chemical potentials and 2dfi their difference. We 
restrict our analysis to zero temperature and work in the mean field approximation. 

The effect of the attractive interaction between fermions is to produce a gap A in the quasiparticle dispersion laws 
which is related to the ipip condensate by 

A(a;) 
< iJaix)'ipp{x) >= -^ e^p . (2) 

The non-homogenous superconductive phase, where the gap is not uniform in space, will be treated in the next Section. 
Here we only discuss the homogeneous case, A(a::) = A = const. In this case the excitation spectrum is described by 
the quasiparticle dispersion laws 

ei = +(5/i + Ep, £2 = -<5/i + Ep , (3) 

with 



Ep = Ve + A2 , ^=py2m-n. (4) 



Using the dispersion laws of the system one evaluates the grand-potential, which is given, at T = 0, by 



^^A2_l f d'p 



G 2 J (2n) 



\ei\ + \e2\-2^ 



(5) 



The integral in this expression is ultraviolet divergent and can be regularized by the usual procedure [3 employing 
the S-wave scattering length a: 

m 1 f (Pp 1 , , 

"" / ^ (6) 



Ana G J (27r)3 p^ ' 

We introduce the dimensionless coupling constant 

nkpa 

where kp is the Fermi momentum. The weak coupling regime, where the BCS approximation holds, corresponds 
to g -^ — oo. This approximation is generally very good for superconductivity in metals. On the other hand, in 
cold atoms the strength of the interaction can be varied working in the vicinity of a Feshbach resonance, where the 
scattering length strongly depends on the applied magnetic field. Therefore both the weak and strong coupling regimes 
can be reached in this case. 

The gap parameter A and the mean chemical potential fi can be determined self-consistently employing the equations 



dil 
dA 

9/i 



= 0, (8) 

(9) 



where n = kp/Sir'^ is the fermionic number density. Eqs. JSJ and (PJ are the gap equation and the number equation 

respectively. In the weak coupling regime the chemical potential fi differs from the Fermi energy ep by an amount of 

order A^//^^; therefore the approximation /i ~ ep is usually adopted and the number equation is not used. 

Let us note explicitly that we do not write equations for fii and fi2 separately. We work at fixed n, but do not 

. . . dn 

impose conditions on 6n = rii — n2. As a consequence, we do not write down a third equation: —-— = —5n, which 

I— 1 . . . ^"^^ . 

would be needed in the analysis if 5n were held fixed |14j . Since Sn is not fixed, conversions between particles of 

different species are allowed and the thermodynamic potential whose minimum defines the vacuum state is not f2, but 

F = n + fm. 

A. Numerical solution of the self-consistent equations 

In order to determine the actual ground state of the system we solve self-consistently Equations ||HJ) and lO to 
determine the values of A and /.( for various values of the coupling g and polarization Sfi. From the computed A and 
fi we evaluate the free-energy F = Q + ^N in the superconductive phase which we compare with the corresponding 
result in the normal unpaired phase. For Sfi = this analysis was performed in Ref. |3ll | and, as a test of our numerical 
code, we have reproduced all the results of this paper. In particular we mention some results valid at Sfi = 0. The 
value of the gap is an increasing function of g. The average chemical potential /i decreases and becomes negative at 
g ~ 0.15 signaling that the system has reached the BEG phase. Finally, the superconductive phase is energetically 
stable. 

This qualitative behavior persists for moderate values of dfi, but increasing Sfi eventually a transition to the normal 
phase occurs. The nature of this transition we now discuss. 

For any fixed value of the coupling constant g we denote by Sfic the largest chemical potential difference that the 
homogeneous superconductive phase can sustain. In other words for Sfi > Sfic the system enters the normal phase. 
Clearly Sfic depends on the coupling strength g. In the weak coupling regime (g — > — cx)) the value of the critical 
polarization approaches the Clogston limit |j| Sfic — Ao/-\/2, where Aq is the value of the gap at Sfi = 0. The phase 
transition between the superconductive and the normal phase is first order. With increasing values of g the value of 
Sfic increases and the phase transition remains first order for values of the coupling smaller of ~ 0.13. 

For values of the coupling in the range 0.13 ^ g ^ 0.175 the gap equation has three non trivial solutions in a range 
of values of Sfi close to the critical mismatch. One of these solutions corresponds to a maximum of the free-energy, 



the other two to local minima. The minima are favored for different values of J/i. At small values of 6^ the favored 
state is the one with A = Aq. For values of Sfj, larger than a critical value the favored state is the second one, with 
A < Aq. The transition between these two states is first order. We remark that such behavior of the free-energy 
takes place only in the range of the coupling 0.13 ^ 5 ^ 0.175. For values smaller than ^ 0.13 there is one phase 
transition from the homogeneous to the normal phase. For values of g large than ~ 0.175 one of the minima of the 
free-energy disappears and, increasing Sfi, one finds a second order phase transition from the normal phase to the 
unpaired phase. 

In order to clarify the behavior in the above-mentioned range of g, we plot in Fig. ^ the free-energy difference 
F ~ Fq {Fq the value at A = 0) as a function of A for various values of Sfi at g — 0.135, i.e. inside the interval 
[0.13,0.175]. For each value of A, the value of /i is determined by the equation dF/dfj, = 0, corresponding to Eq. 
©. We notice that, since the total number density is fixed, the average chemical potentials of the broken (A 7^ 0) 
and normal (A = 0) phases are in general different. For 5fi ~ 0.936 ep the free-energy has a global minimum at 
A = Aq — 0.95 ep and a local minimum at A ~ 0.75 ep] at 6fi — 0.940 ep the two minima are almost degenerate, and 
the values of the gap at the local minima are A = Aq and A ~ 0.625 ep; finally for (5/.t = 0.942 ep the former local 
minimum becomes the global one (and vice-versa), and the gap at the global minimum is A ~ 0.6 ep. For higher 
values of Sfi the value of the gap decreases monotonically and for 6)1 = Sfj^c ~ 0.955 ep the system has a second order 
phase transition to the normal phase. 
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FIG. 1: Free energy difference F — Fq as a function of A for various values of Sfj, at g = 0.135. 



The dependence on g of the order of the phase transitions is shown on the left panel in Fig.[21by three representative 
values of the dimensionless coupling constant, one inside the interval [0.13, 0.175], another one on the left, and a third 
one on the right of the interval. The lowest curve refers to g = —0.1. We have not considered here the possibility of 
inhomogeneous superconductivity and therefore we have a first order phase transition from the superconductive to 
the normal state. It occurs at S^ ~ 0.79Ao. For 0.79 ^ J/i/Ag < 1 the superconductive phase becomes metastable 
and is shown as a dotted line. The highest curve is computed at g = +0.2: for this value the transition from the 
superconductive to the normal phase is second order. The intermediate curve is obtained at g = -1-0.135 and shows, 
in agreement with the results of Fig. ^ two phase transitions: a first order phase transition from the value 0.95 to the 
value 0.65 of the gap parameter, and a second order phase transition to the normal phase. The values corresponding 
to the metastable phases are depicted as dotted curves. An enlarged picture of this case is in the inset. On the right 
panel in Fig. |21we show the behavior of the average chemical potential ^ as a function of (5/i, for the same values of 
the dimensionless coupling constant g. In the figure, the upper curve (green online) represents the average chemical 
potential in the normal phase The inset refers again to g = 0.135. 

It is also worth mentioning that the first order phase transition between the two minima of the free-energy cor- 
responds to a phase transition between a gapped and gaplcss phase. The gapless phase is characterized by having 
one zero in the quasiparticle dispersion law at one sphere in momentum space. Had the dispersion laws two zeros 
then the system could live in the Breached Pairing phase |3, H^ , but this possibility is not realized in this model at 
least within the present approximations. To illustrate this point we have reported in Figure |31 the results for /x/A vs. 
Sfi/A as lines (green online) labeled with various values of g. Since for some values of g there are first order phase 
transitions, some regions of this diagram are never reached by stable physical states, which is why in some cases the 
lines are interrupted. Such regions are above the thick full (red online) line, which has been determined comparing 
the energies of the various phases, and have been labeled with the letters A, B and C. The regions labeled as A and 
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FIG. 2: On the left: The gap A/e^ vs. S^/ep for three values of the dimensionless coupling g. From top to bottom the lines 
refer to t; = 0.2 (purple online), g — 0.135 (red online) and g — —0.1 (black online). The g = 0.2 curve shows a second order 
phase transition to the normal phase at Sn/ep — 1.27. The g = —0.1 curve shows a first-order transition to the normal phase 
at Sfi ~ 0.39eF. The intermediate curve {g = 0.135), shown in more detail in the inset, shows the existence of two phase 
transitions. One phase transition is first-order. It leads to a superconductive phase with a different, smaller, value of the gap. 
The second transition leads smoothly to the normal phase. On the right: t-i/ep vs. Sfi/ep for the same three values of the 
dimensionless coupling g. The continuous upper curve (green online) refers to the normal phase {g — > — cxs). The other three 
curves from bottom to top refer to g = 0.2 (purple online), g = 0.135 (red online) and g = —0.1 (black online). The inset 
represents an enlargement of the curve at (j = 0.135. In both panels the dotted parts of the g = 0.135 and of the g = —0.1 lines 
correspond to metastable states. 

C correspond to metastable points that are local minima of the free-energy. For g ~ 0.135 they were reported in the 
insets of Fig. |21 as dotted points in the upper curve and lower curve respectively. The points in region labeled as B 
correspond to unstable BP points that are maxima of the free-energy. The remaining parts of the diagram correspond 
to allowed regions. The white area corresponds to the stable gapped phase and the shadow area (yellow online) , with 
the exclusion of the region C, to the stable gapless superconductive phase. In the shadow region 5fi > \J y? 4- A^ 
there are gapless excitations at one sphere in momentum space. All the regions meet at the point P, on the line 
corresponding to g = 0.175. The meaning of this point will be clarified below. 
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FIG. 3: Full (green online) lines are simultaneous solutions of the gap and number equations for different values of the coupling 
constant and of the mismatch between the Fermi spheres. The regions above the full (red online) line labeled with A, B and 
C correspond to phases where no physical solutions of the gap and number equations have been found. The phase in region 
B, the Breached Pairing phase, corresponds to unstable solutions. Regions A and C correspond to metastable phases. The 
shadow (yellow online) region, with the exclusion of the region C, corresponds to the region with stable gapless solutions at 
one sphere in momentum space. The remaining part of the diagram corresponds to stable gapped solutions. 



III. LOFF PHASE 

For values of (5/i larger than du^ it can be energetically convenient for the ferniionic system to form Cooper pairs 
with non-zero total momentum Q- In the following we will consider a simple non homogeneous LOFF condensate 
characterized by a single plane wave 



< ipaix)-)pf3{x) >= Aeape 



z2q-x 



(10) 



where 2q is the total momentum of the pair. More complicated patterns, such as those arising by more plane waves, 
may lead to states with a lower free-energy, but for our illustrative purposes the ansatz (|10[l is sufhcient. 
In the LOFF phase the quasiparticle dispersion laws are given by 



ei(q) = +S^i 



where we have defined S,{q) = 



qp 



p + q 
2m 

n{q)^ 



vW 



e2(q) = -5^i + 



q p 

m 



VM' 



— jjL. The free-energy for the LOFF phase can be written as 



G 



<Pp 

(27r)3 



|ei(q)| + |62(q)|-2e(<7) 



(11) 



(12) 



and the integral is regulated as in the homogeneous case, i.e. by employing the S-wave scattering length defined in 
Eq.®. In addition to the gap and number Equations (jH)),© one has to consider the equation 



dq 



0, 



(13) 



which determines the modulus of q; the direction of q is spontaneously determined by the system. Therefore in the 
non-homogeneous phase one has to solve a system of three coupled equations: the number equation, the gap equation 
and Ea. HLJ|) for /x, A and g as a function of g and (5/i and to look for minima of the free-energy F = ^ + /in. 

From general arguments [23 one knows that the one-wave LOFF phase is energetically favored for values of the 
mismatch 5^ in some interval: 5^i < Sfi < 5^2- Let us define the amplitude of the LOFF window as (S/i2i = 6fi2 — SfJ-i. 
Since the free-energy in the one-wave LOFF phase differs only slightly from the value in the normal phase, it turns 
out that 6fii ~ S^,c. In the previous Section we have found that 5fic is an increasing function of g. Comparing the 
free-energy of the LOFF phase with the free-energy of the normal phase we obtain that also dfi2 is an increasing 
function of g but Sfi2i turns out to be a non monotonic function of the coupling. In the weak coupling regime we 
recover the well known result 6112 — 0.754Ao and therefore Sfi2i = (0.754 — 0.707)Ao = 0.047Ao is an increasing 
function of g. However in the intermediate coupling regime 6fii grows faster than Sfi2 and the amplitude of the LOFF 
window, depicted in Fig^ reaches its maximum amplitude at g ~ —0.1 and then begins to shrink; for g ~ 0.05 the 
LOFF phase disappears. 
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FIG. 4: Amplitude of the LOFF window 5fi2i = 5/i2 — 5/ii as a function of tlie coupling g. The LOFF window has its maximum 
amplitude for g ~ —0.1 and shrinks to zero for g > 0.05. 



The largest value of the LOFF window is about 0.035 ej?. However, as we have already mentioned, this value is a 
lower bound for Sfi2i because more complicated crystalline structures might be able to sustain larger values of the 
mismatch between the Fermi surfaces. 



As a check of our results, we have determined the second order phase transition hues of the phase diagram by a 
Ginzburg-Landau (GL) expansion of the grand potential fl, both in the homogeneous and in the LOFF phase. The 
use of this approximation is justified because also in the strong coupling regime one has A/Sfi -^ near the second 
order lines. Since we are interested to the second order transitions, it is enough to expand Q up to the fourth order 
in A, so the grand potential can be written as 



n^no + -a2 + ^A^ , 



a 

— i 
2 



where fJo is the free gas contribution and the coefficients are given by 



n—~oc 



d?p 



(27r)3 (jw„-ei)(iw„ + e2) 



(14) 



(15a) 



d^p 



'-^.RW 



n—'-oo 



(27r)3 {iujn - eiy{iujn + (.lY 



In Eqs. H15|l the e^ are the dispersion laws of the quasi-particles, 



ei 



(p + q) 

2m 



Ml , 



£2 



(p-q) 

2m 



M2 



(15b) 



(16) 



(the homogeneous case is studied by putting q = in the above expressions). The divergence in the integral defining 
the coefficient a is cured, by the introduction of the S-wave scattering length, as discussed in Section^ Using the GL 
expansion we reproduce within a few percent the second order transition lines obtained by the numerical evaluation 
of the free-energy minima in the full theory . 

IV. PHASE DIAGRAM 

We summarize our results in the phase diagram depicted in Fig. [S| In the following discussion of the phase diagram 
we will show that there is a correspondence between some regions and lines of the phase diagram and of the diagram 
depicted in Fig. 13 
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FIG. 5: Phase-diagram at T = 0. The full line (red online) indicates the first order phase transition between the homogeneous 
gapped superfluid phase and the LOFF phase for g < 0.05 or the normal phase for 0.05 ^ g ^ 0.13 or the gapless homogeneous 
superfluid phase for 0.13 ^ 5 ^ 0.175. The dashed line (black online) indicates the second order phase transition between the 
LOFF phase and the normal phase. The dot-dashed line (blue online) indicates the second order phase transition between the 
homogeneous superconductive phase and the normal phase. The dotted (green online) line, which does not correspond to a 
phase transition, separates the homogeneous gapped phase from the homogeneous gapless phase. In the inset it is shown that 
the full line continues beyond the point where the dot-dashed line and the full line meet. 



To begin with we describe the full line (red online). At values of 5 < 0.05 it indicates the first order phase transition 
between the homogeneous superconductive phase and the LOFF phase. For 0.05 ^ g ^ 0.13 the full line indicates a 



first order phase transition from the superconductive phase to the normal phase. The regions corresponding to the 
LOFF phase and to the normal phase are not depicted in Fig. |3| 

For 0.13 ^ g ^ 0.175, the full line represents the first order phase transition between the superfluid gapped phase 
and the superfluid gapless phase (see the discussion in Section ^ in connection with Figs. ^ and [21). From the inset 
of Fig.[Slone can see that such a first order phase transition line terminates at the point P, at g ^ 0.175. This point 
is depicted in Fig. 13 by the same letter. The full (red) line shown Fig. corresponds to the full (red) line shown in 
Fig-El Oiic can understand the different topologies of Figs. [SI and in the following way. Since the regions A, B and 
C in Fig. [21 are not physical they cannot be reported in the phase diagram of Fig. \S\ One can naively say that in the 
phase diagram the areas of regions A, B and C have been shrunk to zero. Therefore once the two branches of the 
full (red) Hue of Fig. [31 are reported in Fig. [31 they overlap as one line. 

For larger values of g ^ 0.13 we do not find a phase transition between the gapped and the gapless phase which 
are separated in Fig. [SI by the dotted (green online) line. Such a line corresponds to the onset of the shadow (yellow 
online) region of Fig. [31 for negative values of /x. The shadow region (with the exception of the region C) shown in 
Fig. [3 corresponds to the region between the dotted green line and the dot-dashed line ( blue onhne) of Fig. [3 The 
dot-dashed line represents the second order phase transition between the superfluid phase and the normal phase. 

The phase transition between the LOFF phase and the normal phase (dashed red line) is of the second order. This 
result, already known in the weak coupling approximation [g, persists also for intermediate values of the dimensionless 
coupling constant, i.e. up to g ^ +0.05. For this value of g and for Sfj. ~ Q.8ep, the transition lines LOFF-BCS 
and LOFF-Normal phase cross. For values of the coupling constant larger than g ^ 0.05 the LOFF phase cannot be 
realized, see the discussion related to Fig. ^ 

In conclusion we have analyzed in the mean field approximation the phase diagram of a fermion superfluid system 
comprising two unbalanced populations. We have worked at fixed total number density but arbitrary density difference, 
conversions between particles of different species are allowed thereof. As a consequence, at T = 0, the phase diagram 
depends on two parameters, the dimensionless coupling constant and the mismatch 6)1 between the two Fermi spheres. 

In view of possible applications to the study of the phase diagrams of systems composed of cold fermion atoms 
the present study should be considered as preliminary. As a matter of fact we have not included two effects. First, 
we have included only fermion self-interactions, without taking into account a confining potential. This development 
is postponed to future investigations. Second, we have not included fluctuations around the mean field solution. In 
the literature this is considered in general as a good approximation, at T = |3ll |. However, as discussed in |32l| . for 
values of the coupling smaller than some critical value (with our definitions, this corresponds to the region \g\ < 0.1) 
one enters a region where fluctuations may play a role. In general including fluctuations has the effect to enlarge 
the region of the ordered phase. This is revealed in the Quantum Monte Carlo computation of [sj by a gap in the 
dispersion law of the quasiparticle in a region where the normal phase should already be present. In other words the 
effect of the fluctuations would be to increase the transition line around g ^ in Fig. [31 Therefore it is plausible that 
a more advanced study, would reveal in the region j^j < 0.1 a more complex structure, more akin to the one proposed 
in [13. We plan to come to this problem in the future. 
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